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I. INTRODUCTION 

With the increasing precision of the measurements of 
the cosmic microwave background (CMB) anisotropics 
and large scale structures (LSS) of the universe as well 
as various other astronomical observations, it is now pos- 
sible to have a clear and consistent picture of the history 
and content of the universe since nucleosynthesis. Al- 
though the now widely used term of "Standard Model of 
Cosmology" might remain premature as our knowledge of 
the cosmological scenario is by far less precise than that 
of the Standard Model of particle physics, the matter con- 
tent of the universe as well as its expansion rate are now 
known within a precision of a few percents with great con- 
fidence. It is also well established that the cosmological 
perturbations which gave rise to the CMB anisotropics 
and the LSS of the universe was inflationary like, with 
a close to scale invariant Harrison-Zeldovitch spectrum. 
Moreover, the measurement of both the temperature and 
polarization anisotropies of the cosmic microwave back- 
ground allows to test the paradigm of adiabaticity of the 
cosmological perturbations and hence the precise nature 
of the mechanism whih has generated them. 

The simplest realizations of the inflationary paradigm 
predict an approximately scale invariant spectrum of adi- 
abatic and Gaussian curvature fluctuations, whose am- 
plitude remains constant outside the horizon, and there- 
fore allows cosmologists to probe the physics of inflation 
through observations of the CMB anisotropies and the 
LSS matter distribution. However, this is certainly not 
the only possibility. Models of inflation with more than 
one field generically predict that, together with this so- 
called adiabatic component, there should also be entropy, 
or isocurvature perturbations 0, 0, 0, 0, IS Eli associ- 



ated with fluctuations in number density between differ- 
ent components of the plasma before photon decoupling, 
with a possible statistical correlation between the adia- 
batic and isocurvature modes 0- Baryon and cold dark 
matter (CDM) isocurvature perturbations were proposed 
long ago [U as an alternative to adiabatic perturbations. 
A few years ago, two other modes, neutrino isocurva- 
ture density and velocity perturbations, have been added 
to the list 0- Moreover, it is well known that entropy 
perturbations seed curvature perturbations outside the 
horizon 0, 0, LI- so that it is possible that a signifi- 
cant component of the observed adiabatic mode could 
be strongly correlated with an isocurvature mode. Such 
models are generically called curvaton models [Tol lllj ]. 
and are now widely studied as an alternative to the stan- 
dard paradigm. Furthermore, isocurvature modes typi- 
cally induce non-Gaussian signatures in the spectrum of 
primordial perturbations. 

In this paper, we describe in more detail the analysis 
performed in Ref [l2T | , constraining the various isocurva- 
ture components. We also extend it by including addi- 
tional observational constraints, and extra free param- 
eters in the model. We use data from the tempera- 
ture power spectrum and temperature-polarization cross- 
correlation measured by the WMAP satellite 01 1 as well 
as from the small-scale temperature anisotropy probed 
by VSA Q, CBI [3 and ACBAR 0; from the matter 
power spectrum measured by the 2-degrcc-Field Galaxy 
Redshift Survey (2dFGRS) [HI and the Sloan Digital 
Sky Survey (SDSS) 0|; and also use data from the re- 
cent type la Supernova compilation of Ref. [l^- We do 
not use the data from Lyman-a forests, since they are 
based on non-linear simulations carried under the as- 
sumption of adiabaticity. The first bounds on isocur- 
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vature perturbations assumed uncorrelatcd modes |2Cj . 
but rece ntly also correlated ones were considered in 
Rcfs. @ El El |H HI 11 Our general analysis 

includes this possibility. In the first part of this work, 
we will not assume any specific model of inflation, nor 
any particular mechanism to generate the perturbations 
(late decays, phase transitions, cosmic defects, etc.), and 
thus will allow all five modes — adiabatic (AD), baryon 
isocurvature (BI), CDM isocurvature (CDI), neutrino 
isocurvature density (NID) and neutrino isocurvature ve- 
locity (NIV) — to be correlated (or not) with each other 
, and to have arbitrary tilts. 

In terms of model building, the simplest situation be- 
yond the paradigm of adiabaticity is that of a single 
isocurvature mode mixed with the adiabatic one. There- 
fore, we shall not consider more than one isocurvature 
mode at a time, and our primordial perturbations will 
be described by three amplitudes and three spectral in- 
dices, associated respectively with the adiabatic, isocur- 
vature and cross-correlated components. This choice is 
somewhat different from that of Refs. 0, 0, |2f|, who 
introduce several modes at a time, but a single tilt for 
each power spectrum of primordial perturbations. The 
assumption that all the modes have comparable ampli- 
tudes and a common tilt are both difficult to motivate 
theoretically and, to our knowledge, all proposed mecha- 
nisms based on inflation stand far from this case. For in- 
stance, double inflation leads to at most one isocurvature 
mode, always with a tilt differing from the adiabatic one; 
on the other hand, the curvaton scenario predicts a single 
tilt, but only one isocurvature mode, fully correlated or 
anti-correlated with the adiabatic mode. However, in the 
absence of any theoretical prior, we believe the approach 
of Refs. 0, 0, |2f| is interesting and complementary to 
ours. 

For simplicity, we will neglect the possible effects of 
spatial curvature and tensor perturbations, and assume 
that neutrinos are massless |2(j. Each model will be de- 
scribed by eleven cosmological parameters: the six usual 
parameters of the standard ACDM model; the amplitude 
and spectral index of the primordial isocurvature pertur- 
bation; the amplitude and spectral index of the cross- 
correlation angle between the adiabatic and isocurvature 
modes; and finally, a parameter w describing the equation 
of state of dark energy, assumed to be time-independent 
as a first approximation. In addition, we will treat con- 
servatively the matter-to-light bias of the 2dF and SDSS 
redshift surveys as two extra free parameters. 

In section II we describe the notation we use in our 
analysis of isocurvature modes and discuss the relation 
with multi-field inflation, and in particular two-field mod- 
els, on which we will concentrate ourselves. In section III 
we discuss the general bounds on our full parameter space 
from CMB, LSS and SN data using a Bayesian likelihood 
analysis. In section IV we analyze those bounds in a con- 
crete model of two-field inflation: double inflation, with 
two uncoupled massive fields. A particular case is that in 
which the two fields have equal masses, like in complex 



field inflation, which we show is not ruled out. In section 
V we present our general conclusions. 



II. NOTATIONS 

For the theoretical analysis, we will use the notation 
and some of the approximations of Refs. 0,0. During in- 
flation more than one scalar field could evolve sufficiently 
slowly that their quantum fluctuations perturb the met- 
ric on scales larger than the Hubble scale during inflation. 
These perturbations will later give rise to one adiabatic 
mode and several isocurvature modes. We will restrict 
ourselves here to the situation where there are only two 
fields, 0i and 02, and thus only one isocurvature and one 
adiabatic mode. Introducing more fields would compli- 
cate the inflationary model and even then, it would be 
rather unlikely that more than one isocurvature mode 
contributes to the observed cosmological perturbations. 

The evolution during inflation will draw a trajectory 
in field space. Perturbations along the trajectory (i.e. 
in the number of e-folds N) will give rise to curvature 
perturbations on comoving hypersurfaces, 
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while perturbations orthogonal to the trajectory will give 
rise to gauge invariant entropy (isocurvature) perturba- 
tions, 



S = 5\n — = 



Spi 



S Pj 
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For instance, entropy perturbations in cold dark mat- 
ter during the radiation era can be computed as iScdi = 
<5cdm — 3(5 7 /4. In order to relate these perturbations dur- 
ing the radiation era with those produced during infla- 
tion, one has to follow the evolution across reheating. 

During inflation we can always perform an instan- 
taneous rotation along the field trajectory and relate 
the gauge-independent perturbations in the fields |27j . 
S<pi = 5<pi + ((j)i/H)ip, with perturbations along and or- 
thogonal to the trajectory, So and Ss, 




cos (p sin ip 
- sin <p cos <p 




(3) 



with ip the rotation angle between the two frames. 

In this case, the curvature and entropy perturbations 
on super-horizon scales can be written, in the slow-roll 
approximation, as 



K k = H 
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where a 2 = <j> 2 + 0|- The problem however is that, con- 
trary to the case of purely adiabatic perturbations, the 
amplitudes of both curvature and entropy perturbations 
do not remain constant on superhorizon scales, but evolve 
with time In particular, due to the conservation of 
the energy momentum tensor, the entropy perturbations 
seed the curvature perturbations, and thus their ampli- 
tude during the radiation dominated era evolves accord- 
ing to iH 



K k =a(t)HS k , 
S k = f3(t)HS k , 



(6) 
(7) 



where a and (3 are time-dependent functions characteriz- 
ing the evolution during inflation and radiation eras. A 
formal solution can be found in terms of a transfer ma- 
trix, relating the amplitude at horizon crossing during 
inflation (*) with that at a later time during radiation, 



(8) 




where the transfer functions arc given by 



Tss{t, t*) = exp 



(3{t')H{t')dt 



(9) 



Tns(t,t*)= / a(t')H(l/)Tss(t',U)dt' . (10) 
Jt, 

Note that in the absence of primordial isocurvature per- 
turbation, £>* = 0, the curvature perturbation remains 
constant and no isocurvature perturbation is generated 
during the evolution. This is the reason for the entries 
T-R.TZ — 1 an d Tsn = 0, respectively, in the transfer ma- 
trix. Note also that in many models, a(t) and (3(i) vanish 
after reheating, so that , ) remain constant during 
radiation domination on super-horizon scales. However, 
this is not true for instance when the fluid carrying the 
isocurvature perturbations has a significant background 
density (compared to the total Universe density), as as- 
sumed in the curvaton scenario. 

Since <pi and <p2 are essentially massless during in- 
flation, we can treat them as free fields whose fluctu- 
ations at horizon crossing have an amplitude 5<pi = 
(-Hfc/V2fc 3 ) ej(k), where H k is the rate of expansion at 
the time the perturbation crossed the horizon (&* = aH), 
and ej(k) are Gaussian random fields with zero mean, 
(e;(k)) = and (e,(k)e*(k')) = %<5(k-k'). Now, 
since 5a k and 5sk and just rotations of the field fluctua- 
tions, they are also Gaussian random fields of amplitude 
Hk/V^k 3 . However, the time evolution (JSJ) will mix those 
modes and will gencrically induce correlations and non 
Gaussianitics. 

Therefore, the two-point correlation function or power 
spectra of both adiabatic and isocurvature perturbations, 
as well as their cross-correlation, can be parametrized 



with three power laws, i.e. three amplitudes and three 
spectral indices, 
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where kg is some pivot scale. Since the time of hori- 
zon crossing in the transfer functions is scale- 
dependent, the correlation angle A(fc) is in general a 
function of k. In the above definitions, we approximated 
cos A(/c) by a power law with amplitude cos Afe and tilt 
n cor . So, we assumed implicitly that the inequality 



( ^- 

k 



< 1 



(12) 



holds over all relevant scales. In the following analysis, 
we will impose that for each value of cos Afc the tilt n cor 
is restricted to the interval in which the inequality holds 
between fc min = 4 x 10 _5 Mpc _1 and fc max = 0.5 Mpc -1 , 
which is roughly the range probed by our CMB and LSS 
data sets. Far from that range, we expect that next- 
order terms become important and that the power-law 
approximation breaks down. The fact that be 
non-zero was already considered in a recent paper |24j . 

The angular power spectrum of temperature and po- 
larization anisotropies seen in the CMB today can be ob- 
tained from the radiation transfer functions for adiabatic 
and isocurvature perturbations, 8f d (fc) and 0} so (fc), 
computed from the initial conditions (TZ Ta d, <S ra( j) = (1, 0) 
and (0,1), respectively, and convolved with the initial 
power spectra, 
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Then, the total angular power spectrum reads 

Ci = A 2 Cf + B 2 Cf ° + 2 A B cos A fco C, cor . (13) 

In many works (see for instance 0,H2), the following 
parametrization is employed: 



Ci = a 
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2 y^iso 
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where the global normalization has been marginalized 
over, and B is the entropy to curvature perturbation ratio 
during the radiation era, B = S ra d/lZ ra d- We will use 
here a sli ghtl y different notation, used before by other 
groups 0,113 , where A 2 = (1 - a) and B 2 = a (up to a 
global normalization factor), so that 

Q = (1 - a) Cf + a Cf° + 2/3 y/a(l-a) Cf OT . (15) 

The parameter a runs from purely adiabatic (a = 0) 
to purely isocurvature (a = 1), while (3 defines the cor- 
relation coefficient, with j3 = +1(— 1) corresponding to 
maximally correlated(anticorrelated) modes. There is an 
obvious relation between both parametrizations: 

a = B 2 /{l + B 2 ), /3 = cosA feo . (16) 

This notation has the advantage that the full parameter 
space of (a, 2(3^/a(\ — a)) is contained within a circle 
of radius 1/2. The North and South rims correspond 
to fully correlated (/3 = +1) and fully anticorrelated 
(/3 = — 1) perturbations, with the equator correspond- 
ing to uncorrelated perturbations ((3 = 0). The East 
and West correspond to purely isocurvature and purely 
adiabatic perturbations, respectively. Any other point 
within the circle is an arbitrary admixture of adiabatic 
and isocurvature modes. 

Note that interpreting the results in non fully corre- 
lated and anti correlated models is made difficult by the 
fact that as soon as n COI ^ 0, the contours one obtains 
depend on the choice of the pivot scale. For example in 
the simple case where n a d = Tii SO , the relative amplitude 
between Cf A and C\ so is unchanged when one changes 
the pivot scale, and hence a is unchanged, but the cor- 
relation angle A(fc), and hence (3 changes with the pivot 
scale as f3(k\) = /3(fco)(fci/fco)™ cor - ^ n our example, if 
n COT > 0, points within the (a, 2/3 -y/ a(l — a)) circle are 
shifted vertically toward the edge of the circle when one 
increases the pivot scale and shifted toward the horizon- 
tal (3 = line when one decreases the pivot scale. 

We should emphasize that the three amplitude param- 
eters ln[10 10 7?. r ad], ot and (3 are defined at k = k , and 
that comparing bounds from various papers is straighfor- 
ward only when the pivot scale is the same. Throughout 
this paper, we will use fcn = 0.05 Mpc -1 , which differs 
from our previous work [l2|. but has the advantage of 
matching most of the literature. This value of k corre- 
sponds roughly to multipole number l Q ~ 300. 

III. GENERAL RESULTS 

We computed the Bayesian likelihood of each cosmo- 
logical parameter with a Monte Carlo Markov Chain 
method, using the public code CosmoMC (28[ with op- 
tion MPI in the COSMOS supercomputer. For each case, 
we ran 32 Markhov chains, and obtained between 30,000 
and 60,000 samples (without including the multiplicity 
of each point) . 



We used the public code CAMB [2^ in order to com- 
pute the theoretical prediction for the Ci coefficients of 
the temperature and polarization power spectra, as well 
as the matter spectra P(k), for all four different com- 
ponents (The BI mode is just a rescaling of the CDI 
mode, see below). Moreover, we slightly modified the 
interface between CAMB and CosmoMC in order to in- 
clude the cross-correlated power spectra, as well as three 
independent tilts — in the CosmoMC jargon, the three 
tilts and the three amplitudes were implemented into the 
code as "fast parameters" , in order to save a considerable 
amount of time. The likelihood of each model was then 
computed using the detailed information provided on the 
experimental websites (or directly from the correspond- 
ing routines of the CosmoMC code when available), using 
three groups of data sets: i) CMB data: 1398 points from 
WMAP (2 < I < 900) [13, 8 points from VSA (580 < 
I < 1700) [3, 8 points from CBI (700 < I < 1800) 0, 7 
points from ACBAR (920 < I < 1960) 0; ii) LSS data: 
32 points from 2dFGRS (up to fc max = 0.1 h Mpc" 1 ) £3, 
17 points from SDSS (up to k max = 0.15 h Mpc" 1 ) pj : 
and iii) SNIa data: 157 points from Riess et al. [l!j . We 
have checked that the inclusion of the HST |3l| and BBN 
priors are irrelevant for the determination of the Hubble 
parameter and the baryon density; that is, the data from 
CMB, LSS and SNIa is enough in order to determine 
these parameters, and therefore we ignored the priors in 
the final analysis. 

In order to take into account the /3-dependent con- 
straint on ricor, see eq. (|12fl . we choose to define a new 
parameter 

tfcar^ncor/Wr 1 ). ( 17 ) 

whose boundaries are fixed once and for all by the values 
of the pivot scale and the scales (fc m i n , fc max ) introduced 
in the previous section. The basis of cosmological pa- 
rameters used by the Markov Chain algorithms consists 
on: (i) the seven parameters describing the standard adi- 
abatic ACDM model, extended to dark energy with a 
constant equation of state, and (ii) the four parameters 
describing the admixture of one (correlated) isocurvature 
mode, already defined in the previous section. More ex- 
plicitely, we use the following basis: 1) the overall nor- 
malisation parameter ln[10 10 7?. r ad] where 7?. r ad is the cur- 
vature perturbation in the radiation era; 2) the adiabatic 
tilt n a d; 3) the baryon density Ub = ^Bh 2 ; 4) the cold 
dark matter density LJ c dm = ^cdmh 2 ; 5) the ratio 9 of 
the sound horizon to the angular diameter distance mul- 
tiplied by 100; 6) the optical depth to rcionization r; 
7) the dark energy equation of state parameter w; then, 
for the isocurvature sector; 8) the isocurvature contri- 
bution a; 9) the cross-correlation parameter (3; 10) the 
isocurvature tilt n- lso ] 11) the parameter <5 cor which deter- 
mines the cross-correlation tilt n cor ; and finally, two ar- 
bitrary bias parameters associated to the 2dF and SDSS 
power spectrum. Our full parameter space is therefore 
1 3-dimcnsional . 

We did not devote a specific analysis to the case of the 
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FIG. 1: The one-dimensional likelihood functions for our basis of eleven independent cosmological parameters (not including 
the tilts of the two redshift surveys), for the adiabatic mode alone (AD) or mixed with the three different types of isocurvature 
modes (AD+CDI, AD+NID, AD+NIV). The first seven parameters are those of the standard ACDM model, extended to dark 
energy with a constant equation of state. The last four parameters (a, j3, rii so , 5 COT ) describe the isocurvature initial conditions. 
(<5 C or is define in Eq. (El)) 
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FIG. 2: Continuation of figure showing the ID likelihood of some derived cosmological parameters, for the same cases. 
These likelihoods should be considered with care, because the parameters shown here do not belong to the basis used by the 
Markov chain algorithm. Therefore, the shape of the above likelihoods depends not only on the likelihood of the underlying 
parameters, but also on the properties of the functions relating them to the parameters of the basis. This explains for instance 
why n C or = S COI ln(|/3| _1 ) seems to be well-constrained, while S COT and (3 are not. 



baryon isocurvature modes, which is qualitatively simi- 
lar to that of CDI modes, since the spectra are simply 
rcscalcd by a factor ^e/^cdm (^B/^cdm f° r the cross- 
correlation): thus, compared to the AD + CDI case, sig- 
nificantly larger values of a will be allowed in the AD 
+ BI case. Like in other recent analyses, we find that 
the inclusion of isocurvature modes does not improve sig- 
nificantly the goodness-of-fit of the cosmological model, 
since in the AD+CDI, AD+NID and AD+NIV cases the 
minimum % 2 is always between 1672 and 1674 for 1614 
degrees of freedom, to be compared with 1674 for 1618 
degrees of freedom in the pure adiabatic case. There- 
fore, the question is just to study how much departure 
from the standard picture is allowed, by computing the 
Bayesian confidence limit on the isocurvature parame- 
ters. A more detailed analysis of model comparison with 
Bayesian Information Criteria will be done in a follow 
up paper. 

On Fig. we plot the marginalised likelihood for our 
basis of eleven cosmological parameters, in the cases 
AD+CDI, AD+NID and AD+NIV, compared with the 
pure adiabatic case. Figure |21 shows the likelihood of 
some derived parameters. It appears that most parame- 
ters are robust against the inclusion of isocurvature per- 
turbations (this is in agreement with the conclusion of 
Ref. |25j| that with only one isocurvature mode present, 
no significant parameter degeneracy pops out). Our 95% 
C.L. on a in the three cases is given in Table [|] 

Note that in the limit a = 0, the three parameters j3, 
fiiso, S cot become irrelevant. So, the fact that pure adi- 
abatic models are very good fits implies that these pa- 
rameters are loosely constrained. This explains why the 



TABLE I: The one-dimensional 2-a ranges on the isocurvature 
mode coefficients for the various models. 



model 


a 2/3 [a(l -a)] 1/2 


AD+CDI 
AD+NID 
AD+NIV 


< 0.6 -0.7 to 0.3 

< 0.4 -0.2 to 0.8 

< 0.3 -0.4 to 0.6 



corresponding likelihoods on Fig. are not wcll-pcakcd 
like for other parameters. In addition, these likelihoods 
should be considered with great care, because it is dif- 
ficult for the Markhov Chains to explore in detail the 
tails of the multi-dimensionnal likelihood corresponding 
to tiny values of a, where basically any value of (j3, rii so , 
5 COT ) are allowed. Therefore, increasing the number of 
samples would tend to flatten these likelihoods, while 
the other ones would remain stable (as we checked cx- 
plicitely). However, it is clear that all models prefer a 
large isocurvature tilt and saturate the bound nj SO < 3 
that we fixed in the present analysis. This feature is 
important for understanding our results and comparing 
with other analyses, as explained in the last paragraph 
of this section. 

On Fig. |3| we plot the two-dimensional confidence lev- 
els directly for the isocurvature and cross-correlation co- 
efficients (a, 2f3y/a(l ~ a)) in the three cases AD+CDI, 
AD+NID and AD+NIV. The last plot corresponds to 
the AD+CDI case with a prior n cor = which is rele- 
vant for the bounds on double inflation, but the results 
are not substantially different from the general AD+CDI 
case. We see that the AD+CDI model slightly prefers 
anti-correlated cases (note that (3 < means a positive 
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FIG. 3: The 2-cr contours of a and the cross-correlated mode coefficient 2f3^J a(\ — a), for a) the CDI isocurvature mode; b) 
the NID mode; c) the NIV mode; d) the CDI mode, with the constraint n CO r = and the contours of equal 2(R 2 — l)/sk = 
0, ±0.5, ±1, ±2 from double inflation. 



contribution from the cross-correlated component), while 
AD+NID and AD+NIV models clearly prefer correlated 
ones. 

On Figure we plot the CMB and LSS power spec- 
tra for two particular CDI and NID models. In order 
to get a better understanding of our bounds, we chose 
models with large values of a, still allowed at the 2-cr 
level: respectively, a = 0.53 and a = 0.41. The detailed 
values of other cosmological parameters for these models 
are given in the figure caption. In the CDI example, one 
can see that the non-adiabatic contributions to all spec- 
tra remain tiny, excepted for the matter power spectrum 
on scales k > 0.2/iMpc -1 , due to the large isocurva- 
ture tilt rij so = 2.93 of the model. This is an indication 
that our a bound in the CDI case depends very much on 
constraints on the small-scale matter power spectrum, 
while future improvements in the determination of CMB 
spectra would not reduce it dramatically. Note that a 
precise experimental determination of the amplitude pa- 
rameter er 8 (which is mainly sensitive to scales around 
k = 0.2 h Mpc -1 ) would probably not change things ei- 
ther, since our CDI models have roughly the same as val- 



ues as purely adiabatic well- fitting models (see Figure|5J); 
on the other hand, any constraint on smaller wavelengths 
could improve the bounds. This means in particular that 
our choice not to include the Lyman-a data plays a cru- 
cial role in our results. 

The same conclusions apply to the NID model of Fig.0] 
In addition, in the NID case, we see that the non- 
adiabatic contribution is significant also for small-scale 
(/ > 200) temperature and polarization spectra. Indeed, 
it is well-known that the NID isocurvature and cross- 
correlated modes can mimick the adiabatic CMB spectra 
to a better extent than CDI modes (essentially because 
the amplitude of the secondary peaks is not strongly su- 
pressed). So, in the NID case, future improvement in 
the CMB data should help to improve the bounds on the 
isocurvature fraction. 

Our bounds are difficult to compare with those of 
Ref. j2^| , because of our different parameter space, ob- 
servational data set and conventions of normalization. In 
the AD + CDI case, our analysis is closer to the one of the 
WMAP team > although we have more free cosmolog- 
ical parameters (n cot , w), more data (SDSS, CBI, VSA) 



and less constraints on the matter-to-light bias. The 
95% limit B < 0.33 obtained by WMAP would corre- 
spond to a < 0.1 in our notations, which is significantly 
smaller than our results. Also, The WMAP la bounds 
on niso are 1.26 ± 0.5, while we find that the likelihood 
peaks at our maximum allowed value rn so = 3. The most 
likely explanation is that the use of the Lyman-a data in 
the WMAP analysis eliminates all our well-fitting mod- 
els with riiso > 2 and large a values. Similar conclusions 
apply to our previous results [TJ, in which we did not 
use any Lyman-a information, but adopted a flat prior 
0.6 < 7ii SO < 1.5 (this was the interval in which our grid 
of models was computed). Then, most of the well- fitting 
models of ^3] had slighly negative values of n a d — ni SD . 
Therefore, translating our previous results in terms of a 
pivot scale ko — 0.05 Mpc -1 would lead to a small de- 
crease in the a bounds, making them comparable with 
the WMAP bound in the CDI case, and smaller than the 
conservative bounds of this work. 



IV. DOUBLE INFLATION 

In this section we will use the previous bounds to con- 
strain the parameters of a concrete model of inflation 
called double inflation 0, Hj|. This is an inflationary 
model with two massive fields, <p\ and (f>2, of different 
masses, with ratio R = mi/m2- No coupling (except 
gravitational) between fields is assumed. The equations 
of motion and Friedmann equation can be written as 



in +ZH4>i+m{<t)i =0, 

l 2 + 3H <j> 2 + mjfo =0, 



(18) 
(19) 



As inflation proceeds, s decreases and 9 also decreases. 
We will call sh = 60 to the number of e-folds before 
the end of inflation when the scale corresponding to our 
Hubble radius today exited during inflation. 



A. Linear perturbations 

The great advantage of double inflation is that it is 
possible to find explicit formulas for the perturbations 
on superhorizon scales. The growing mode solutions for 
the scalar fields and the scalar metric perturbation in the 
longitudinal gauge, and in the slow-roll approximation, 
are given by 



-Ci(*0 



H 

IP 



C 3 (k) 



2{R 2 - l)R 

3(R 2 cb 2 + ( b 2 ) 2 



6<fn = Ci(fc) 

H2 = <Mfc) 

60 H 



C 3 (k) 

cm 



2H<% 



R 2 4>\ 
2HR 



2^2 



R 2 <b 2 



(23) 
(24) 
(25) 



Since <j>i and <f>2 are independent uncoupled scalar fields 
and essentially masslcss during inflation, we can use the 
general formalism of section II, and write, in the slow-roll 
approximation, 



Ci(fc) 



C 3 (k) = - 



H, 



3 (0!ei(k)+0 2 e 2 (k) 



3H k /ci(k) e 2 (k) 



2V2fc 3 \m\(bi m\<b2 



(26) 



(27) 



H 2 
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<Pi + 4>l + m 2 cb 2 +m 2 cb 2 ), (20) 



where k 2 = 8irG. In the slow-roll approximation, <b 2 <C 
mfeb 2 ; fa <C Hfa, we have the solutions j3^| 



-2) **i 



4.s 



1/2 



cos ( 



with s = — ln(a/a en d) the number of e-folds to the end 
of inflation. For simplicity, we restrict ourselves to the 
case where (b\ and </> 2 remain positive during inflation. 
Substituting into the rate of expansion, 

H 2 {s) ~ ^m 2 2 s (l + (R 2 - 1) sin 2 0] , (21) 

and, using the equations of motion, one can integrate out 



where Hk is the rate of expansion when the perturbation 
of wavenumbcr 



k H e s 



(28) 



left the horizon during inflation, where the scale of our 
present horizon is = 3000 ft- -1 Mpc. 

We will now assume that the light scalar field decays 
at the end of inflation into the ordinary particles, giving 
rise to photons, neutrinos, electrons and baryons, while 
the cold dark matter (CDM) arises from the decay of the 
heavy field. In principle, part of the CDM could also be 
produced by the light field or the heavy field could also 
decay into ordinary particles, but we will ignore these 
possibilities here. Then, the perturbations in the comov- 
ing gauge take the form 



<W 
ds 



00 



R 2 



1 



tan 6* 



2s 1 + R 2 tan 2 9 



(22) 



n B 



n B 



to obtain 



and there is only one isocurvature mode, the CDI mode, 
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FIG. 4: Temperature, E-polarization and matter power spectra for two particular CDI and NID models. In order to get a 
better understanding of our bounds, we chose here two models with large values of a, still allowed at the 2-a level: respectively, 
a = 0.53 and a = 0.41. Other parameter values are for the CDI (resp. NID) model: uj b = 0.0217 (0.0196), w cdm = 0.112 
(0.131), 9 = 1.06 (1.01), r = 0.068 (0.131), w = -0.88 (-1.44), n ad = 0.96 (1.02), n iso = 2.93 (2.95), n cor = 0.05 (0.03), 
ln[10 10 ft rad ] = 3.73 (3.98), (3 = -0.62 (0.88). The \ 2 of the two models is respectively 1675 and 1674. From top to bottom, 
we show the Cj T , C™ and P(k) power spectra, as well as the contribution of each component: adiabatic, isocurvature, 
cross-correlated, total. The CDI isocurvature and cross-correlated components have been rescaled by a factor indicated in each 
figure. We also show the data points that we use throughout the analysis, from WMAP (black), ACBAR (grey), CBI (blue), 
VSA (yellow), 2dF (black) and SDSS (blue). In the case of the matter power spectrum, one should not trust a "x 2 -by-eye" 
comparison with the data: first, because the spectrum has to be convolved with the experimental window function before the 
comparison (this changes its slope significantly); second, because we show here the data points before rescaling by the two bias 
factors, which are left arbitrary for each experiment. 
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all of which are gauge invariant quantities. During the 
radiation era, the initial conditions of all these modes are 
described in terms of only two ^-dependent quantites, 
and Sk- The pure adiabatic initial conditions are given 
by the gravitational potential during the radiation era, 

$rad(fc) = ^rad(fc) = |ci(fc), 

with C\{k) the amplitude of the growing adiabatic mode 
during inflation. On the other hand, the isocurvaturc 
initial conditions in the radiation era arise from the per- 
turbations in the heavy field 4>\ at the end of inflation. In 
the long wavelength limit, the perturbations of this field 
during reheating follow closely the field itself, so that 
its energy density perturbations satisfy, in the comoving 
gauge, 



Sp 



(c) 



p 



which is constant during inflation, across reheating and 
into the radiation era. The entropy perturbation is dom- 
inated by the CDM density perturbation during the ra- 
diation era, <S ra d ~ ^dm- Using the values of C\{k) and 
Cs(k) during inflation, we can finally write 



"fcrad(fc) 



Srad(k) 




sin0 fe ei(k)+cos0 fc e 2 (k)) ,(29) 
ei(k) i? 2 e 2 (kr 



sin tffe 



cos 9 k 



(30) 



Note that the cross-correlation amplitude cos A(fc) is al- 
most independent of k, due to the cancellation of the 
factor Sfe between 7?. ra d(fc) and 5 ra d(fc), and to the fact 
that 9k is a mild function of k, see Eq. (I22|) . In this case, 
the tilt n cor vanishes. The adiabatic and isocurvaturc 
tilts can be computed explicitly using (|28|l 
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1 



(1 + (R 2 + R 4 ) tan 2 9 + R 6 tan 4 i 
^From Vv, and Vs, we obtain the parameters 
R 4 tan 2 9 + 1 



s 2 k sin 2 







(R 2 



R 4 tan 2 9 
1) sin 6» 



1 



Vi? 4 tan 2 9 + 1 



(32) 

(33) 
(34) 
-(35) 

(36) 
(37) 



We have plotted these parameters as a function of the 
angle 9 in Fig. 03 for s/. ~ 60. 




FIG. 5: The values of parameters a and (3 as a function of 
the angle 8, for different values of the ratio R = rai/m2 in 
double inflation. 

These parameters have maximum and minimum values 
given by 



(38) 
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as can be seen in Fig. [S] 



B. Bounds on double inflation 

In this subsection we will impose the general bounds 
found in section III to the double inflation model |33j . 
Note that in this model it is assumed that the heavy field 
decays into cold dark matter, and therefore we only have 
one isocurvature component, CDI. The bounds from CDI 
will be used to constrain this particular model. We will 
leave for the future a detailed analysis of other two-field 
models of inflation. 

In order to derive specific constraints, it will be useful 
to take into account the following relation between a and 
(3, see Eqs. JSgJ and P7|) 



2/V"(l-") = 



R 2 



1 



Sk/2 



(1-a), 



(40) 



which corresponds to a straight line in the contour plot 
of Fig. 3. This way, one can evaluate the likelihood at 
which a given value of R is ruled out. Unfortunately, 
from the contours in Fig. 3, one cannot restrict much the 
range of R, except to exclude R > 5 at 2er. 

Even if the model passes this constraint for a given 
_R, it is possible that the prediction on n ac i (|33|l and nj so 
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(|35|l do not agree with the bounds on these parameters. 
In our case, the bounds are so loose that any tilt is al- 
lowed. Perhaps in the future, with better observational 
constraints, we may use the information on the tilts to 
further rule out double inflation models. 



C. Massive complex field 

Another model worth exploring is the particular one 
in which the two fields have equal masses, corresponding 
to a massive complex field $ = -^("fii + i4>i)' If we 
rewrite $ = -^a exp(i(f), with modulus a and phase ip, 
the lagrangian can be written as 

Note that there is no potential for the phase, so it will be 
free to fluctuate, which will induce a large isocurvature 
component, as we will see, and which can be used to rule 
this model out. 

In this case, the curvature and entropy perturbations 
are 



"KradO) = ^Lv'sfc e CT (k) , 




e s (k). 



V2fc 3 y/Sk 

with e CT and e s orthonormal. Therefore 
1 



4 + 1 



p = o. 



(41) 
(42) 

(43) 



The curvature perturbation has a tilt n at \ = 1 — 2/s = 

0. 97, but the isocurvature perturbation has no tilt, n; so = 

1, and the two modes are uncorrelated, (3 = 0. For the 
moment, this model is not ruled out. 



the inclusion of an isocurvature component, see Table 1 
and Fig. 3. So, the pure adiabatic scenario remains the 
most economic and attractive scenario. 

In contrast with the WMAP analysis, we decided not 
to include any data from Lyman-a forests, since con- 
straints on the linear power spectrum coming from these 
experiments are derived under the assumption of a plain 
adiabatic ACDM scenario. We did not include either 
strong priors on the isocurvature spectral index, unlike in 
our previous work, and allowed this parameter to vary up 
to rii so = 3. This conservative approach leads to a prefer- 
ence for models with a very blue isocurvature primordial 
spectrum, and to upper bounds on the isocurvature frac- 
tion significantly larger than in other recent analyses: on 
a pivot scale k = 0.05 /iMpc -1 the amplitude of the cor- 
related isocurvature component can be as large as about 
60% for the cold dark matter mode, 40% for the neutrino 
density mode, and 30% for the neutrino velocity mode, 
at 2(7. This leaves quite a lot of freedom, for instance, 
for double inflation models with two uncoupled massive 
fields. Assuming that one of these fields decay into Cold 
Dark Matter, our results simply imply that the mass of 
the heavy field cannot exceed five times that of the light 
field at the 2a confidence level. It is expected l35l that 
in the near future, with better data from Planck and 
other CMB experiments, we will be able to reduce further 
a possible isocurvature fraction, or perhaps even discover 
it. The present results also suggest that constraining the 
linear matter power spectrum on scales which are mildly 
non-linear today will also be crucial in this respect. 



Note added 



After completin g th is work, we came upon the paper 
of Parkinson et al. |22|, where a similar analysis was done 
for a class of hybrid inflation models. 



V. CONCLUSIONS 

Using the recent measurements of temperature and po- 
larization anisotropies in the CMB by WMAP, together 
with recent data from VSA, CBI, and ACBAR; the mat- 
ter power spectra from the 2dF galaxy redshift survey 
and the Sloan Digital Sky Survey, as well as the recent 
supernovae data from the SN Search Team, one can ob- 
tain stringent bounds on the various possible isocurva- 
ture components in the primordial spectrum of density 
and velocity fluctuations. We have considered correlated 
adiabatic and isocurvature modes, and find no significant 
improvement in the likelihood of a cosmological model by 
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